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1.0 SUMMARY 


An important question facing the designers of efficient, high performance aircraft is whether the poten- 
tial advantage of laminar flow can be realized. Advances in materials technology, most notable the ma- 
turing of composite applications, indicate that the surface smoothness requirements for laminar flow are 
now attainable. However, there are only limited data available on the effects of wing sweep, Reynolds 
number, Mach number, and pressure distribution shape on transition behavior and location. 

The objective of this experiment was to develop and use hot film devices to measure transition on a 
Cessna Citation III business jet at various flight conditions and correlate the results with linear bound- 
ary layer stability analyses. These results would provide valuable additional data in establishing criteria 
for use in designing high-performance swept wings that promote long runs of laminar flow. 

This document describes the method and gives the results of the boundary layer stability analyses of 18 
conditions flown during the flight tests. The transition point varied from 5% to 35% chord for these condi- 
tions, and both upper and lower wing surfaces were included. Altitude varied from 10,000 to 43,000 ft and 
Mach number from 0.3 to 0.8. Four cases were at nonzero sideslip. 

Although there is much scatter in the results, the analyses of boundary layer stability at the 18 condi- 
tions led to some of the same conclusions as those from the earlier F-lll laminar flow glove flight tests. 
The three upper surface cases that involved sideslip showed that crossflow and Tollmien-Schlichting dis- 
turbances may interact and reduce the maximum amplification of either isolated disturbance that can be 
tolerated without causing transition. However, the majority of the results indicated that crossflow insta- 
bilities were the primary cause of transition in these tests; the crossflow amplification factors were scat- 
tered around 5 at transition while Tollmien-Schlichting amplification factors varied from 0 to 8. The 
lower surface shows much lower Tollmien-Schlichting amplification at transition than the upper surface, 
but similar crossflow amplifications. No relationship between Mach number and disturbance amplifica- 
tion at transition could be found. The quality of these results is uncertain because of questionable wing 
surface quality, inadequate density of transition sensors on the wing upper surface, and an unresolved 
pressure shift in the wing pressure data. 

The results of this study show the need for careful preparation for transition experiments. Preparation 
should include flow analyses of the test surface, boundary layer disturbance amplification analyses, and 
assurance of adequate surface quality in the test area. The placement of necessary instruments and use- 
fulness of the resulting data could largely be determined during the pretest phase. 



2,0 INTRODUCTION 


2.1 BACKGROUND AND OBJECTIVES 

An important question facing the designers of efficient, high-performance aircraft is whether the poten- 
tial advantage of laminar flow can be realized. Advances in materials technology, most notable the ma- 
turing of composite applications, indicate that the surface smoothness requirements for laminar flow are 
now attainable. However, there are only limited data available on the effects of wing sweep, Reynolds 
number, Mach number, and pressure distribution shape on transition behavior and location. Reference 1 
reports results of a study of sweep effect on transition location for a range of Reynolds numbers, but at a 
Mach number less than 0.3. Transition locations from several general aviation aircraft flight tests are 
reported in reference 2 and include a good range of Mach, Reynolds number, and wing sweep. Neither of 
these references involves a correlation of actual transition locations obtained with the anticipated 
growth of boundary layer instabilities. This connection was made in a NASA EET sponsored flight test of 
a natural laminar flow glove on the F-lll TACT aircraft. Linear stability analyses were correlated with 
transition locations inferred from boundary layer measurements in an attempt to establish transition 
criteria which could be used for Natural Laminar Flow (NLF) wing design. Reference 3 explains the 
results and methods used in that program. However, the determination of transition location by measur- 
ing the boundary layer profile after transition introduces an uncertainty into the work of reference 3. 
Recently, experimentation with hot film devices has provided a new tool for sensing transition. 

The objective of this experiment was to develop and use hot film devices to measure transition on a high- 
performance business jet at various flight conditions and correlate the results with linear stability analy- 
ses such as those used in reference 3. These results would provide valuable additional data in establish- 
ing criteria for use in designing high-performance swept wings which promote long runs of laminar flow. 


2.2 APPROACH 

The aircraft used in this experiment was the Cessna Citation III (Model 650). This aircraft has a quarter 
chord sweep of 25 deg and is capable of operation up to Mach 0.8 and altitudes over 45,000 ft. A wing area 
5 ft in span and 60% in chord was smoothed with fiberglass and filler to eliminate the roughness of fas- 
teners and joints. No attempt was made to alter the shape of the wing sections. 

To determine the pressure distribution on the test area, flights were made with three pressure belts: one 
on each side of the test area and one in the center. On later flights, the center belt was removed and re- 
placed with a series of hot film transition-sensing gages. The correlation of the pressure distributions 
from the outer belts with the middle belt distribution was used to infer the pressure distribution over the 
transition sensors for the second series of flights. The flight test conditions ranged from 10,000 to 43,000 
ft and from Mach 0.3 to 0.8, and positive and negative sideslip was also flown. 

The flight conditions and wing pressure data were recorded on Cessna Aircraft digital flight-recording 
equipment. The hot film signals were also recorded on the same magnetic tape on analog tracks. Signals 
from several of the hot films could be monitored simultaneously in flight with oscillographs. Data reduc- 
tion and hot film signal processing occurred at a ground station. 
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Examination of the hot film data led to the selection of conditions which received detailed analysis. The 
transition location and corresponding wing pressure distributions for the selected flight condition were 
used to determine the boundary layer instability growth up to the transition point. A finite-element infi- 
nite swept wing boundary layer analysis code was used in this process. These boundary layer parameters 
were then used in a linear stability analysis code, which calculated the growth of both Tollmien- 
Schlichting and crossflow disturbance modes. Plots of the disturbance growth that occurred before transi- 
tion give an indication of the amount of growth which could occur on other wings (new designs) before 
transition occurred (i.e., a transition criterion). 
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3.0 SYMBOLS AND ABBREVIATIONS 


Acronyms 


CF 

crossflow 

NLF 

natural laminar flow 

TS 

Tbllmien-Schlichting 

WBL 

wing butt line 

Symbols 



A 

A„ 

c 

C P 

c f 

Cl 

C n 

s 

H 

h 

M 

N 

Re! 

Re e a .i. 

rj.e. 

5 
s 

W 

X 

z 

a* r 

r s 

6 
A 
A 


disturbance amplitude 

disturbance amplitude at neutral stability point 
chord 

wing streamwise pressure coefficient 
surface friction coefficient 
section lift coefficient 
section normal force coefficient 
section pressure coefficient 
altitude 

height of surface above mean profile 
Mach number 

disturbance amplification factor (In A/A 0 ) 
unit Reynolds number 

attachment line Reynolds number based on boundary layer momentum thickness 

leading-edge radius of a wing section 

estimate of standard deviation 

distance measured along a surface 

aircraft weight 

distance measured along the chord line 

distance measured perpendicular to the chord line 

component of dimensional wavenumber in the direction of wing sweep 

boundary layer thickness 

wing sweep angle 

wavelength 
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T ell thickness ratio of the ellipse that best fits the wing section shape at the leading edge 

y disturbance propagation angle, relative to local potential flow direction, positive toward 

wing outboard 

co disturbance frequency 

Subscripts 


a 

c 

n 

T 

oo 


average 

chordwise 

normal to the sweep 
at boundary layer transition 
freestream conditions 
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4.0 FLIGHT TEST PRESSURE DATA ANALYSES 


4.1 METHOD OF ANALYSIS 

The boundary layer stability analysis of the flight test data in the present study requires a detailed 
knowledge of the boundary layer parameters on the wing section of interest, especially near the leading 
edge. Since these could not be measured practically, they were calculated by a boundary layer code. The 
boundary layer code used a finite difference formulation to calculate laminar or turbulent boundary lay- 
ers on two-dimensional, axisymmetric, or infinite swept bodies. As input, the program requires certain 
freestream and geometry parameters and the pressure distribution on the body or wing. Since a founda- 
tion of the stability analysis is the pressure distribution, an accurate and detailed knowledge of the pres- 
sure distribution was important in this study. 

As mentioned previously, the first series of test flights measured pressure distributions at three sections 
in the test area: WBL 111, 141, and 171. The primary section of interest, 141, was bracketed by the other 
two so that when the middle pressure belt was removed for the transition-sensing flights, the two outside 
belts could give an indication of the pressure at the primary section. At each flight condition of interest, 
several complete pressure distributions for each section were taken. These data were reduced to pressure 
coefficients, and all distributions for each flight condition and wing section were plotted together. Any 
obviously incorrect data points were apparent on these plots; after eliminating them, the remaining 
points from each wing pressure orifice were averaged. This resulted in one pressure distribution per 
flight condition for each wing section, as well as an indication of the variation in the distribution that 
might occur while trying to hold the flight condition. The relationship between the pressure distributions 
of the two outside sections (WBL 111 and 171) and that of the primary section (WBL 141) was one of the 
guides in determining the pressure distributions for the primary section during the later transition- 
sensing flights. The pressure measured on the primary section during the pressure flights could not be 
used directly on the later transition-sensing flights because the flight weights were not the same and the 
accuracy and detail of the measured distribution was not adequate. 

A general layout of the instrumentation arrangement for this study is shown in figure 1. Three 
scanivalves were used to measure wing section pressures, and a detailed description of the components 
and their function is in reference 4. 

Another guide to determining the pressure distributions of the primary section was two-dimensional 
transonic flow analyses of the primary section’s airfoil shape. The forward 60% of this airfoil was deter- 
mined by careful measurements taken after the flight test. Since transition occurred near the leading 
edge for almost all flight conditions, it was important to have a fine resolution of the pressures near the 
leading edge. The computer analysis was a better method of finding these details than interpolating from 
the measured wing pressures on the sides of the test area. 

Using these theoretical results and the measured pressure distributions from the outside sections during 
the transition-sensing flights led to the pressure distributions to be used for the boundary layer analyses. 
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4.2 PRESSURE DATA RESULTS 


The method of taking wing pressure data during this study resulted in a large amount of data. As many 
as 25 individual pressure distributions for each of the three wing sections at one flight condition were 
obtained. Averaging these distributions was appropriate in order to reduce the volume of data and be- 
cause of the small degree of scatter in the pressures at each flight condition. Figure 2 shows a set of pres- 
sure distribution data consisting of 10 samplings for one of the flight conditions and each of the three 
wing sections. These data were taken in about 45 sec (10 pressure readings per second per scanivalve). 
The small amount of scatter is readily apparent. Estimates of standard deviation of pressure coefficient, 
S Cp , at a given test condition, which were calculated during the pressure averaging process, were gener- 
ally less than 0.005. In areas where the pressure was very sensitive to flight condition, such as near the 
wing leading edge and near shocks, the S Cp rose to about 0.02. The few cases of interruption of the flow of 
data, due to reasons which were not determined, were easily spotted by looking at the estimates of stan- 
dard deviation. Those few points were then averaged manually. There were a few pressure orifices on the 
wing belts which gave incorrect readings (fig.2), and these were omitted from further processing. Figure 3 
shows the averaged pressure coefficient data for wing sections 111, 141, and 171 from the pressure- 
measuring flights. These graphs were used as a guide to establish the difference in pressure distributions 
between the test section (WBL 141) and the outside sections. 

The three-dimensional boundary layer and stability codes are set up to operate on a normal (perpendicu- 
lar to the sweep) section of the wing. Therefore, the streamwise pressure data from the transition-sensing 
flights were transformed into normal pressure coefficients by dividing the averaged C p ’s by the cosine 
squared of the local wing sweep angle (simple sweep theory). Figures 4 to 21 show the resulting aver- 
aged, then “normalized”, pressure coefficient data from wing sections 111 and 171 for all cases analyzed 
during the transition-sensing flights. Table 1 outlines the streamwise conditions for each of the 18 cases 
selected for stability analysis. 

A two-dimensional transonic viscous airfoil analysis code (ref. 5) also played an important part in deter- 
mining the pressure distribution to be used in the boundary layer stability analyses. The airfoil coordi- 
nates to use with the two-dimensional analyses came from two sources. Over the front 60% of the chord, a 
detailed and accurate measurement of the wing at WBL 141 was made after the flight tests. Aft of this 
area the coordinates of the wing design loft were used initially. These two sets of coordinates were com- 
bined smoothly, and this airfoil was then scaled up by 1/cos A to account for the streamwise to normal 
transformation (simple sweep theory). The sweep used here and later for the stability analyses was 26.6 
deg, which is the wing sweep at about 7.5% chord. This value was chosen because most of the transition 
points during this study were between 5% and 15% chord. This scaled up airfoil could now be analyzed by 
the two-dimensional code using normal Mach number and lift coefficient for any of the cases and the 
results compared to the “normalized” flight data for that case. This procedure showed that the airfoil 
first used for the analyses had too much aft camber, resulting in too little lift being carried forward. Sev- 
eral variations in the airfoil shape were tried, some of which had extended chords to simulate a separa- 
tion bubble, since wing pressures near the trailing edge are low, indicating separation. The one that was 
most successful and used for the remainder of the study was the first airfoil used with a small amount of 
the lower surface cusp filled in, but no extended chord. Figure 22 shows a trace of this airfoil. Even after 
the airfoil was found which gave a chordwise lift distribution that compared well to the data, there was a 
pressure coefficient level shift between the theoretical analysis and data. This shift was similar to that 
which would arise from using an incorrect static pressure measurement during data reduction, but a 
careful review of the data acquisition and reduction procedure did not locate any errors. This shift in 
pressure coefficients was also evident when comparing the flight data with unpublished wind tunnel data 
from NASA Langley Research Center 8-ft Pressure Tunnel Test 795. The closest comparison possible 
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with the wind tunnel data was pressures from WBL 160 at Moo = 0.70. Figures 3m and 3n both show the 
comparison, since Moo = 0.70 falls between the Mach 0.68 and 0.725 flown during this program and 
shown on those figures. The sectional lift coefficient of the flight and wind tunnel data were within a few 
percent. The wind tunnel configuration did include the nacelles,' so the shift was probably not caused by 
the presence of the nacelles or engine power setting, and an estimate of the body pressure field indicates a 
much smaller effect than the observed discrepancy. The c p level shift was also not likely to be caused by 
translating the measured wing pressure coefficients to equivalent normal section pressure coefficients, 
which comes only from simple sweep theory. This procedure is used often at Boeing on wings of this class 
of sweep and taper with good results. Since no definite cause for the Cp discrepancy could be found, the 
theoretical pressures were shifted to best match the data for each of the 18 cases analyzed. This is illus- 
trated in figure 4. Thus, the theory was only really used for interpolating the experimental pressure data. 

To summarize, the “normal” section conditions used to calculate the theoretical pressure distributions to 
compare to the data were defined as: 


• M n = Moo cos (26.6) 

° c i n ~ c n /cos 2 (26.6) 

© Reynolds number was calculated based on M n , test altitude, and a normal chord length of 5.42 ft. 
Table 2 gives the normal section conditions for the 18 study cases. 

Figures 4 to 21 also show the theoretical pressure distributions plotted with the corresponding normal 
section data for all 18 cases. The value of the shift necessary to match the two is reported for each. The 
irregularities in the theoretical distributions near the leading edge also appear in the curvature distribu- 
tion of the WBL 141 airfoil definition. The apparent discrepancies between the data and theory arise be- 
cause of the influence of the pressure flight data. Reviewing figure 3 points out that the middle section of 
the test area appears to have a slightly different section shape, and the greater density of pressure ori- 
fices also gives greater insight into the details of the rapid pressure decrease near the leading edge, 
which is a characteristic of this airfoil. The theoretical pressure distributions as shown in figures 4 to 21 
are the distributions that were used in the infinite swept wing boundary layer code for the stability anal- 
yses. 

Figure 23 summarizes the data analysis procedure from flight pressure data through boundary layer sta- 
bility analysis. 
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5.0 BOUNDARY LAYER STABILITY ANALYSES 


5.1 Stability Analysis Method 

The method selected for analyzing boundary layer transition data in this study is based on linear bound- 
ary layer stability theory (ref.6 to 13). The basic premise of this method is that transition is caused by the 
amplification of initially small disturbances in the boundary layer. When the amplitude of the distur- 
bance becomes sufficiently large, nonlinear effects cause transition to turbulent flow. Although the later 
stages of transition are beyond the scope of the theory, it still provides the best currently available basis 
for correlating transition data. 

The use of stability theory for correlating transition data was first used in detailed calculations by Smith 
(ref. 14), whcwfound that the total theoretical amplification from the neutral point to the transition point 
was approximately e 9 for a large number of two-dimensional data. In principle, this discovery allowed 
transition to be predicted for any airfoil by computing the theoretical amplifications. This has been gen- 
erally confirmed by later research, but the exponent has been found to range from 8 to 15. For swept 
wings, however, transition can occur much sooner than predicted by the e 9 method. It is now known that 
there are several types of instability that must be considered. 

On a high-speed swept wing there are four basic types of laminar boundary layer instabilities to be con- 
sidered (refs.l and 15): (1) Tollmien-Schlichting (TS), (2) crossflow (CF), (3) Taylor-Goertler(TG), and (4) 
leading-edge attachment line stability. 

Tbllmien-Schlichting instability relates to the velocity profile in the local freestream direction. Amplifi- 
cation of TS disturbances is usually small or negative in regions of favorable pressure gradient and large 
in regions of adverse pressure gradient. For incompressible flow, the most unstable disturbance waves 
propagate in the freestream direction, but for compressible flows the waves may be oblique. 

Crossflow instability refers to waves propagating in a direction nearly perpendicular to the local frees- 
tream direction. On swept wings, a crossflow develops within the boundary layer so that the profile of a 
velocity component in directions oblique to the local freestream has an inflection point. Such velocity pro- 
files are highly unstable, so that crossflow instabilities may be dominant on swept wings. Boundary layer 
crossflow is most severe in the wing leading-edge and trailing-edge regions, where pressure gradients are 
largest. 

Taylor-Goertler instability occurs in the flow over concave surfaces. Because the Cessna natural laminar 
flow (NLF) glove did not have concave surfaces in the region designed to have laminar flow, this type of 
instability was not considered in this study. 

Attachment line instability refers to the behavior of the boundary layer along the forward stagnation or 
attachment line (i.e., the locus of points for which the chordwise velocity is zero). The boundary layer flow 
along the attachment line can be either laminar or turbulent depending on Reynolds number and envi- 
ronment, as described in reference 15. If the attachment line flow does become turbulent, the flow over 
the wing will be turbulent also. Detailed stability calculations are not used to assess attachment line 
instability; previous experimental studies have shown that if the attachment line momentum thickness 
Reynolds number is less than 100, the leading-edge flow will be laminar. In this study, the attachment 
line Reynolds number was somewhat below 100 at the transition test section. 
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The empirical determination of the allowable theoretical amplification factors for swept wings in flight is 
an objective of the present analysis. The rate by which a disturbance is amplified depends not only on the 
local flow conditions and boundary layer profile, but also on the disturbance frequency and propagation 
direction, so there are several possible variations of the process. The process used here is outlined below, 
and reference 3 contains details of the method as actually applied. 

In situations in which both TS and CF instabilities are present, as in the present case, interactions may 
occur that cause earlier transition than if either was absent. Such interaction was evident in the F-lll 
laminar flow flight test, reference 3. Following Smith (ref. 14), it is more common to work with the natu- 
ral logarithm of the amplification factors rather than their actual values. The logs of the amplification 
factors are usually called N-factors. For the F-lll data, amplification factors for both modes were calcu- 
lated, and each experimental transition point plotted on a the N TS - N CF plane as shown in figure 24. The 
shaded band representing the data gives an indication of the effect of such interactions. The data band of 
figure 24 can be used as a transition criterion. The transition data from the present study were analyzed 
in a similar manner to assess and improve the F-lll criterion. 

The boundary layer temperature and velocity profiles, which are the primary output of the boundary 
layer code, became the primary input to the stability program, which is a Boeing modification of a compu- 
ter program originated by L. M. Mack (ref. 16). This program solves the boundary layer stability equa- 
tions for three-dimensional, linearized, parallel flow for a perfect gas and can calculate either Tollmien- 
Schlichting or crossflow stability. The program was used to calculate disturbance growth curves, such as 
those shown in figure 24. 

The disturbance growth direction used to compute amplification factors was along the local potential flow 
streamline. Mack (ref. 11) determined that this was a satisfactory approximation to the actual growth 
direction, which is equal to the real part of the group velocity angle. The stability code can compute ei- 
ther incompressible or compressible stability. In the present study, compressible stability was used for all 
but one of the calculations. The sixth-order equations (which neglect dissipation) were used instead of the 
complete eighth-order equations. This use resulted in a significant reduction in computation time, and, 
as shown by Mack (ref. 11), results given by the sixth-order equations for a transonic swept wing differed 
from those of the eighth-order equations by only a few percentage points. 

TS disturbances were followed downstream, keeping frequency and wave angle fixed. The wave angles 
changed from one case to the next and were determined by examining the change in amplification for 
different wave angles for each case. The angle was then chosen for each case to correspond closely to the 
angle for maximum disturbance amplification. For each case several frequencies were calculated to es- 
tablish an envelope of TS amplification. 

CF disturbances were followed downstream, keeping the frequency fixed and letting the wave angle vary 
in accordance with the irrotationality condition applied to the wavenumber vector, as proposed by Mack 
(ref. 11). The result of applying the irrotationality condition to an infinite swept wing analysis is that the 
spanwise (in the direction of the wing leading edge) component of the dimensional wavenumber, a* rs , 
must remain constant as the disturbance propagates downstream. Therefore, in defining the envelope of 
CF disturbances, disturbances having a range of a* rg values are followed downstream. The frequency was 
kept at zero for this study, although zero frequency waves may not be the most highly amplified. This 
method of analysis is suggested by flow visualization of swept wing transition, which often shows 
streamwise striations indicating the presence of uniformly spaced streamwise vortices. The analysis just 
described corresponds to this flow model. 
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5.2 STABILITY ANALYSIS RESULTS 


Eighteen cases (Table 1) were chosen for stability analysis from the large volume of flight data recorded 
during the transition flights. Most of the cases are from data taken at 30,000-ft altitude, cases 7 to 18, 
with four cases from 10,000 ft, cases 1 to 4, and only two from the high altitude of about 43,000 ft, cases 5 
and 6. All cases were taken with the aircraft at near zero yaw except the last four cases. Cases 15 and 16 
were taken at a negative yaw angle (nose right) and cases 17 and 18 at a positive yaw. Normal chord 
Reynolds number for the 18 cases ranged from a low of 5.42 x 10 6 at the high altitude, case 6, to a high of 
13.8 x 10 6 at the low altitude, case 4. Normal Mach number ranged from 0.264 for cases 1 and 2 to 0.715 
for case 14. Thble 2 summarizes the normal flow conditions and calculated instability growth parameters, 
N ts and N cf , near the transition point for each of these cases. 

The boundary layer analysis of case 1 indicated separation (c f = 0) at a surface distance aft of the attach- 
ment line, s/c, of 0.0725. This is not surprising when one considers the strong pressure peak near the 
leading edge for this condition (fig. 4). The measured transition location for this case was at or ahead 
of x/c = 0.05, which is very close to an s/c of 0.09 for this condition. Therefore, it is likely that a laminar 
separation bubble caused the transition in this instance. Linear stability calculations would be of no use 
in predicting this mode of transition because the extremely adverse pressure gradient causes almost in- 
stantaneous growth in TS instability. Therefore, none were made for case 1. It is important to note here 
the significance of the difference between chordwise distance from the leading edge, x/c, and surface dis- 
tance aft of the attachment line, s/c. The surface distance was used in the boundary layer analyses and 
subsequent stability analyses, whereas the transition gage locations were measured chordwise aft of the 
wing leading edge. Because the instability growth was rapid for many of the cases analyzed, it was impor- 
tant to know the transition gage locations with respect to the attachment line location for each case. Fig- 
ure 25 shows a plot of s/c versus x/c with attachment line locations for the 18 cases. 

For all cases except case 1, the boundary layer code did not predict laminar separation for at least 10% 
chord beyond the point where the transition was measured. The stability calculations were therefore 
meaningful for these cases and were carried out, when possible, to a point 5% to 10% chord aft of the point 
where transition was measured. For several cases, the crossflow instability calculations were difficult 
and could be completed only to about the point of transition. This occurred on upper surface cases where 
the pressure peak was followed by a rather strong adverse pressure gradient; figure 6, case 3, is an exam- 
ple. Figures 26 to 42 show instability growth plots for cases 2 through 18. Both Tbllmien-Schlichting (TS) 
and crossflow (CF) disturbance growth factors (N) are plotted for each case in which there was growth. 
The strong favorable pressure gradient for the lower surface, cases 2 and 7, led to no TS growth before 
transition, and the other lower surface cases, 5, 6, 9, 13, 16, and 18, showed very little TS amplification. 
The TS amplification plots show the growth of disturbances of various frequencies, as well as the enve- 
lope of the growth. The wave angle used for the TS calculations for each case is noted. This angle was the 
one which gave the greatest N factors and was determined for each case by calculating the growth for 
different angles. For crossflow disturbance growth, the range of the value of spanwise wave numbers 
used in the calculation is noted on the N CF - s/c plots. For crossflow zero frequency was always used 
(standing waves only), as described above. 

The lower surface cases had strong favorable pressure gradients for the first 20% to 30% chord, so TS 
disturbance growth was insignificant except for cases 5 and 6, which had laminar runs of about 30%. For 
these lower surface cases, the CF disturbances were the ones to be more highly amplified, precisely be- 
cause of the strong favorable pressure gradients in combination with the high wing sweep. The upper 
surface cases showed strong TS growth starting at about 5% to 6% chord. This was the area where the 
strong favorable gradient at the nose ended in a peak minimum pressure followed by an adverse gradi- 
ent. For these cases the CF disturbances showed most of their growth in the first 5% to 8% chord and 
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were generally damped aft of there. Figures 43 and 44 show N TS versus N CF traces for cases 2 to 18 as the 
boundary layer grows aft of the attachment line. The dramatic growth in TS disturbances aft of s/c = 0.06 
for the upper surface cases is easily seen in these plots. 

The hot film transition sensors were located at 2.5% and 5% chord and then at every 5% chord back to 
70% on the upper surface and 65% on the lower surface. From intermittancy data, which is described in 
reference 4, transition occurred rapidly, with each gage generally indicating either fully laminar 
(~ 15% intermittancy) or fully turbulent (~ 85% intermittancy) flow. On the traces of N TS - N CF , the 
location in s/c of the last gage that indicated laminar flow and the first gage that indicated turbulent flow 
can be noted. It can be said that transition occurs at a combination of TS and CF disturbance growths 
somewhere between these two points. These two points and the trace between them are plotted for each 
case in figures 45 and 46. Unfortunately, in general, there was a large difference in N-factors between 
the last laminar point and the first turbulent point, resulting in a poorly defined transition band. This 
was basically because transition generally occurred far forward where the N-factors grew rapidly in the 
short distances between the sensors. 

Figure 45 is quite different in character from the F-lll data (fig. 24). One interpretation of this figure 
would be that crossflow instability was the primary cause of transition, with TS having little influence. 
For this highly swept wing, this may be the case. Figure 46 also supports that conclusion, but in this case 
there was almost no TS amplification. The three sideslip cases, 6, 8, and 10, did show some interaction 
(fig. 47) but much less than the F-lll data. The equivalent lower surface cases, 7, 9, and 11, still show no 
TS growth and little change of N CF at transition for the different sweep angles. The low TS N-factors are 
the natural result of the lower surface pressure distribution and not inconsistent with the upper surface 
data. Transition generally occurs at N CF ’s of 4 to 6 for both surfaces. 

Apparently low N TS factors at transition for the lower surface were also noticed in the F-lll NLF glove 
flight test (ref. 3) (fig. 48). Again, the low TS values are due to the pressure distribution, but the real 
question is why the low N CF at transition. The reason for these low values on the lower surface was not 
determined for that test, although noise contamination from the F-lll underwing inlet was a possibility. 

Since transition generally occurred earlier than anticipated on the basis of the F-lll data, explanations 
were sought. Potential culprits are waviness and attachment line contamination, or unknown Mach 
number effects. Figure 49 shows a waviness trace taken from measurements of the wing profile at WBL 
141. This trace was calculated by taking the perpendicular distance between the measured coordinates 
and a smoothed set of the same measured coordinates. The large spike near the leading edge is artificial 
and due to an error in reference points used during the wing measurement. The other waves shown in 
figure 49 may or may not present a problem for sustaining laminar flow. From reference 17 an acceptable 
wave criterion of 


2h ^ T 59 000 (cos 2 A)T/2 

T < L A Re, 1 5 J 


is established for single waves and one-third of this for multiple waves. Note the h as defined herein is 
half the amplitude as used in reference 17. For the most stringent flight condition (lowest altitude, high- 
est speed), case 4, the allowable wave height, h, for multiple waves is about 0.0017 in for a 2-in wave- 
length or 0.0025 for a 5-in wavelength. The test surface meets this requirement only very marginally and 
definitely does not at a few sharp spikes like those at which s equals 50, 53, and 69 in; therefore, surface 
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waviness cannot be ruled out as a possible adverse influence on laminar flow during this study. If these 
waves produced initial disturbances in the boundary layer, perhaps the lower speed flow on the lower 
surface amplified them more, accounting for the low N TS factors at transition on that surface. 

As for attachment line contamination, the criterion is summarized in reference 15 as: 


Re 


^ 0.405 sin A 

>/ COS A’ 




(1 + ^ 11 .) 


(r ' e j2 < 90 


For this study, at the highest Re loo (low altitude, high speed), case 4, the Re 0al is about 90 according to 
this formula. These calculations indicate little possibility of a problem due to attachment line contamina- 
tion and transition at the section where transition was measured. If the attachment line boundary layer 
becomes turbulent ahead of the line of transition gages, one would expect that all gages would register 
turbulence. However, if disturbances arose near, but on either side of the attachment line inboard of the 
transition gauge location, they could be convected outboard as they moved downstream. This would be 
more pronounced on the lower surface where the chordwise component of velocity is lower. 

The data in this study were taken from Mach numbers of 0.3 to 0.8, but as shown in figure 50, the N- 
factors at transition from the boundary layer analyses show no obvious dependence on Mach number. The 
N cf data show a definite grouping around 4 to 6, independent of Mach number, but the N xs data show no 
discernable grouping. 
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6.0 CONCLUSIONS AND RECOMMENDATIONS 


6.1 CONCLUSIONS 

The analysis of the boundary layers and their stability for the 18 cases supplied to Boeing from the Cita- 
tion III laminar flow flight test resulted in the following conclusions. 

• The quality of the transition data from the hot film sensors during this study seems good. However, 
there are uncertainties which cloud the usefulness of the stability analyses. First is the lack of mea- 
sured pressure data on the wing along the section where transition was measured, leading to heavy 
reliance on theoretically determined pressure distributions for the stability analyses. The shift that 
was necessary to match analytical and experimental wing section pressures was also never ex- 
plained and may point to an error in the experimental pressures. Second, the rapid growth of TS 
disturbances on the upper surface made the 5% chord spacing of the hot films almost too sparse. 
Third, the quality of the surface in the area where transition was measured was marginal with re- 
gard to waviness and was not adequately surveyed. Although the attachment line Reynolds number 
was probably not too high, it would have been desirable to isolate the inboard leading edge from the 
leading edge in the test area. 

With these reservations, the following results were obtained: 

• Transition generally occurred sooner than expected on the basis of the reference 3 data. 

• There is little evidence of an interaction between the TS and crossflow instabilities. Since such an 
interaction is to be expected, the three sideslip cases were examined separately. These data do show 
an interaction, but much weaker than for the F-lll, and most of the cases indicate that crossflow 
instability apparently caused transition, independently of the magnitude of Tollmien-Schlichting 
instability growth. 

• A wide range of Mach numbers was flown during these tests and no dependence on freestream Mach 
number could be determined for TS or CF N-factors at transition. 


6.2 RECOMMENDATIONS 

The results of this study point strongly toward the need for detailed preparation before conducting 
laminar flow experiments: 

• The pressure distributions that will be encountered during the experiment should be determined 
and verified beforehand. This will influence: 

• whether the experiment will be useful without revisions. 

• whether a change in the wing section shape would make the experiment much more useful 
(decreasing the “peakiness” of the section in this study would have been very helpful). 

• the density and placement of instrumentation. 
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• The quality of the surface throughout the area where laminar flow may occur needs to be very high 
and verified before laminar flow testing begins. If practical, a good way to do this would be to use 
flow visualization at one or two of the more ambitious (high transition Reynolds number) test condi- 
tions. 

• It would be very helpful to have pressure and transition measurements on the same surface simul- 
taneously. 

• Suitable density of transition sensors is necessary if the results are to be correlated with theoretical 
stability analyses. 

• If a laminar flow experiment may involve high sweep angles, the span of the test area should be 
considerable, and the attachment line boundary layer inboard of the test area should be diverted so 
there is a fresh boundary layer on the attachment line in the test area. 

• The results of this experiment call into question the existence of a well-defined transition boundary 
on the Nrpg vs N CF diagram. Additional analyses should be made of other possible correlations or 
other methods of applying the stability theory. 
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Table 1 - Flight Conditions For Cases Analyzed 


Case 

Flight 
. Number 


Altitude, 

ft 

Weight, 

lb 

Static 

Temperature, 

op 

Surface 

xlc (First 
Turbulent 
Gage) 

Sideslip, 

deg 

1 

1180 

0.295 

9965 

16765 

20 

Upper 

0.05 

0 

2 

1180 

0.295 

9965 

16765 

20 

Lower 

0.10 

0 

3 

1180 

0.491 

9980 

16576 

19 

Upper 

0.10 

0 

4 

1180 

0.538 

9976 

16060 

18 

Upper 

0.05 

0 

5 

1182 

0.442 

29918 

15133 

— 

Lower 

0.15 

0 

6 

1182 

0.632 

29837 

15411 

-53 

Upper 

0.10 

0 

7 

1182 

0.632 

29848 

15377 

-53 

Lower 

0.10 

0 

8 

1182 

0.631 

29884 

15034 

-53 

Upper 

0.15 

-3.2 

9 

1182 

0.635 

29848 

14981 

-52 

Lower 

0.15 

-2.5 

10 

1182 

0.635 

29931 

15004 

-53 

Upper 

0.05 

4.3 

11 

1182 

0.631 

29874 

14962 

-52 

Lower 

0.05 

3.9 

12 

1182 

0.679 

29797 

16093 

-52 

Upper 

0.15 

0 

13 

1182 

0.724 

29782 

15523 

-53 

Upper 

0.10 

0 

14 

1182 

0.771 

29679 

15863 

-52 

Upper 

0.05 

0 

15 

1182 

0.770 

29692 

15816 

-52 

Lower 

0.05 

0 

16 

1182 

0.799 

29645 

15703 

-53 

Upper 

0.05 

0 

17 

1181 

0.581 

42696 

15871 

-56 

Lower 

0.35 

0 

18 

1181 

0.680 

42716 

15589 

-59 

Lower 

0.30 

0 




Table 2 - Stability Analysis Conditions and Results 






Cj (Used in 
Airfoil 

Analysis Code) 

A (Used in 
Stability 
Analysis), deg 

At First 
Turbulent Gage 

Case 

Surface 

M n 

Re c x 10 6 

N ts 

N C f 

i 

Upper 

E2 1 

7.60 

0.850 

26.6 

— 

— 

2 

Lower 

m 

7.60 

0.850 

26.6 

0 

5.1 

3 

Upper 

0.439 

12.60 

0.406 

26.6 

8.65 

6.0 

4 

Upper 

0.482 

13.80 

0.338 

26.6 

1.75 

6.05 

5 

Lower 

0.396 

6.10 

0.840 

26.6 

0 

6.05 

6 

Upper 


8.67 

0.475 

26.6 

7.15 

3.15 

7 

Lower 


8.67 

0.475 

26.6 

0 

4.6 

8 

Upper 

i 

8.89 

0.448 

23.4 

8.05 

2.3 

9 

Lower 


8.89 

0.448 

24.1 

0 

4.35 

10 

Upper 


8.40 

0.514 

30.9 

2.25 

5.4 

11 

Lower 


8.40 

0.450 

30.5 

0 

4.4 

12 

Upper 

0.601 

10.90 

0.450 

26.6 

7.4 

3.8 

13 

Upper 

0.648 

9.92 

0.388 

26.6 

3.75 

5.35 

14 

Upper 

0.689 

10.57 

0.350 

26.6 

0.25 

4.65 

15 

Lower 

0.689 

10.57 

0.350 

26.6 

0.4 

5.2 

16 

Upper 

0.715 

10.80 

0.353 

26.6 

0.2 

4.9 

17 

Lower 

0.520 

4.82 

0.944 

26.6 

0.8 

7.85 


18 Lower 0.609 5.40 0.722 26.6 0.5 6.05 
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Figure 1 Instrumentation Layout 
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Fig. 3c Citation III Averaged Wing Pressures - Pressure Flights (continued) 
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Fig. 3n Citation III Averaged Wing Pressures - Pressure Flights (continued) 
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FIG. 23 Overall Stability Analysis Procedure 
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Fig. 25 Relationship of Surface and Chordwise 
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Fig. 26 Citation III Stability Analysis Results, Case 2, Lower Surface 
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Fig. 27 Citation III Stability Analysis Results, Case 3, Upper Surface 
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Fig. 29 Citation III Stability Analysis Results, Case 5, Lower Surface 
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Fig. 30 Citation III Stability Analysis Results, Case 6, Upper Surface 
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Fig. 31 Citation III Stability Analysis Results, Case 7, Lower Surface 
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Fig. 32 Citation III Stability Analysis Results, Case 8, Upper Surface 
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Fig. 33 Citation III Stability Analysis Results, Case 9 Lower Surface 
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Fig. 34 Citation III Stability Analysis Results, Case 10, Upper Surface 
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Fig. 35 Citation III Stability Analysis Results, Case 1 1, Lower Surface 
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Fig. 36 Citation III Stability Analysis Results, Case 12, Upper Surface 
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Fig. 38 Citation III Stability Analysis Results, Case 1^ ,LJpper Surface 
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Fig. 39 Citation III Stability Analysis Results, Case 15, Lower Surface 
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Fig. 40 Citation III Stability Analysis Results, Case 16, Upper Surface 
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Fig. 47 Effect of Sideslip on Instability Growth at Transition 
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